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FOREWORD 


This report was prepared by the Research and Engineering 
Division, Aerospace Group, of The Boeing Company under 
Air Force Contract F33615-72-C-1580, "Propulsion System 
Installation Corrections", Project 1366. The program was 
concucted under the direction of the Prototype Division, 
Air Force Flight Dynamics Laboratory, Air Force Systems 
Command. Mr. Gordon Tamplin was the Air Force Program 
Monitor. 


The program was initiated on 31 December 1971 and draft 
copies of the final reports were submitted for approval 
on 31 October 1972. 


Dr. P. A. Ross was Program Manager and Mr. W. H. Ball 

was principal investigator during the program. Significant 
contributions to the program were made by the following 
individuals: Mr. Joe Zeeben, engine performance; Dr. Franklin 
Marshall, iniet and exhaust system technolegy; Mr. John Petit, 
nozzle internal performance and nozzle/afterbody drag; and 

Mr. Gary Shurtleff, progrzmming. 


This report contains no classified information extracted 
from other ‘\assified reports. 


Publicatic: <£ this report does not constitute Air Force 
approval . «e report's findings or conclusions. It is 
publishe} , for the exchange and stimulation of ideas. 
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Chief, Prototype Division 
Aix Force Flight Dynamics Laboratory © 
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ABSTRACT 


This report presents the results of a research program to 
develop a procedure for calculating propulsion system in- 
Stallation losses. These losses ijiaclude inlet and neazle 
internal losses and external drag losses for a wide variety 
€ of subsonic and supersonic aircraft configurations up to 
Mach 4.5. The calculation procedure, which was Largely 
developed from existing engineering procedures and axperi-+ 
Mental data, is suitable for preliminary stucies of 
advanced aircraft confiaurations. Hndgineerina descriptions, 
equations, and flow charts are provided to help in adapting 
& the calculation procedures to digital commuter roatines. 
Many of the calculation procedures have already bean bregr om ced 
on the CNC 6600 computer. Program listings and flow chart. 
are provided for the calculation preeedures that. have been 
programmed. The work accomp) ished during the wroyy ae ds Cae 
| tained in four separate volumes. Volume I contains an 
¢ engineering description of the calculation procedures. 
Volume II is a programmer's manual containing flow cnarts, 
listings, and subroutine descriptions. Volume IL! contains 
Sample calculations and sample inpnt data. Volus TY con- 
tains fF ookkeeping definitions and data correlations. 
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SYMBOLS AND NOMENCLATURE 


D Drag, 1b.; Hydraulic Diameter, a » in., diameter, in. 


¢ 
A Area, in” 
A, Inlet capture area, in? 
C Ay Local stream tube area ahead of the inlet, in? 
Ay Free-stream tube area of air’ tering the inlet, in? 
Bg 
R Aspect ratio, dimensionless 
« 
B Velocity, decay exponent, dimensionless 
Cc Sonic velocity, ft/sec. 
C,, Drag coefficient, D » dimensionless 
¢ dAper 
C=-D Convergent-divergent 
Cc Additive drag coefficient,'C = Dapp » dim nsionless 
Papp Dapp GA. 
.¢ . c 
Cy Afterbody drag coefficient, ie » dimensiorvless 
A X 
— (PLP) Apage 
¢ Cy Base drag coefficient —————--——=—' climensionlzss 
Base X : 
c Scrubbing 4 fficient, DRAG.» gimensionie- 
Ds ing drag coe Client, : imensionle:s 
¢ F 
Cy Thrust coefficient, ——L— , dimensionless 
g w 
——(V 
g ( cp! 
Cy Stream thrust coefficient, dimensionlcss, | 
( (defined by Figure 48 of Volume IV) | 
Cy Nozzle velocity coefficient, dimensionless | 
Conv. Convergent | 
{ 
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SYMBOLS AND NOMENCLATURE (Continued) 


Thrust, lb. 


Ideal gross thrust (fully expanded), ib. 


Fuel/air ratio, dimensionless 


Gravitational constant, ft/sec 


Enthalpy per unit mass, BTU/lb.; height, in. 
Enthalpy of fan discharge flow, BTU/lb 


Enthalpy of primary exhaust flow after heat addition 
BTU/1b 


Throat height, in? 


Velocity decay coeffucient, dimensionless 


Length, in. 


Mach number, dimensionless 


Pressure, lb/in? 


Relative pressure; the ratio of the pressures P, 
and P, corresponding to the tomperaturcs T. awa 
respectively, along a given isentrope, aieednloiice: 


Total pressure, lb/in? 


Effective heating value of fuel, BTU/lb. 


Dynamic pressure, 
Radius, in. 
Nozzle centerline 
Temperature, °R 


Velocity, ft/sec 


lb/in? 


spacing, in. 
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SUBSCRIPTS 


B 
b, base 


BP 
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SYMBOLS AND NOMENCLATURE (Continued) 


Mass flow, lb/sec 


‘Bleed air removed from engine, 1b/sec. 


Corrected airflow, lb/sec. 


Weight flow rate of fuel, lb/sec. 


Weight flow rate of air, primary plus 
secondary, lb/sec. 


Primary nozzle airflow rate, lb/sec. 


Length, in. 
Angle of attack: convergence angle of nozzle, degrees 


Ratio of specific heats, dimensionless 


Diffuser ioss coefficient, 


, dimwnensionless 
Burner efficiency, dimensionless 


Kinematic viscosity, ft?/sec. 


Density, lb/ft’ 


Burner 

Base flow region 
Bypass 

Boundary layer bleed 
Boattail 


Core (nozzle); capture (inlet) 


Design conditions 


SYMBOLS AND NOMENCLATURE (Continued) 


SUBCRIPTS ae) 
e Boattail trailing edge 
: EFF Based on effective area 

ENG Refers to engine demand 3 

j f£ Fuel 

g gross 

: GEOM Based on geometric area 5 
int Interference; internal 
ip Ideal, primary exhaust flow 
jet Exhaust flow jet 
max Maximum 
min Minimum 
8 Scrubbing flow region ‘ 
SPILL Spillage 
T Total condition; throat station 
ty Total condition, fan flow ‘ 
1/0 Takeoff | 
ty Total condition, primary flow 
° Local conditions ahead of inlet 
2 Compressor face station 
8 Nozzle throat station 
9 Nozzle exit station 
18 Fan discharge throat station 
Free-stream condition 
x Local 
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SECTION I 


GENERAL DESCRIPTION OF COMPUTER PROGRAM 


The engine installation calculation procedure naturally 
breaks into three divisions: inlet, engine, and nozzle/ 
afterbody. Each of these divisions is represented in 

TEM-333 by a set of subprograms logically subordinate to a 
monitoring subprogram which functions as a quasi-main 

program for that set. The term module will be used to 
designate a quasi-main and its subordinate programs. Further, 
within a module procedures are subprogrammed to isolate the 
procedures from one another and to isolate logical dir:ction 
pi 29gramming from computational programming. The benefits 

of this organizetion are threefold. First, checkout. Lach 
module can be checked out independently of the others. Secondly, 
logical simplicity. The logical structure of the program 

as a whole is more apparent. Individual subprograms are 
simpler; being small with a limited purpose. Finally, ease 

of modification. Engine installation programs are constantly 
undergoing modification so the ability to attach and detach 
new modules or procedures quickly is a necessity. The small 
subprogram-module arrangement allows for quick change since 
the connection to the next level is more apparent. 


The modules are tied together by subroutine CQMPUTE which 
calculates the installed performance for a single uninstalled 
performance point. The main program, TEM-333, handles calls 
to the input routines, then reads the uninstalled performance 
points one=by~one; calling C@MPUTE for each, outputting data 
when required until the last point has been computed, 
Programs included in TEM=-333 and their purposes are: 

TEM=333 - Main Program, monitors program logic 

COMPUTE - Controls calculation of a performance point 

Inlet Module 

TABIN = Reads Inlet Input Deck 

INLDRAG ~ Monitors Inlet Drag Computation 


AIRBYP ~ Computes External-Compression Mode 
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AIRSPL - Computes Mixed Compression Mode 
SIZINL - Computes Inlet Capture Area 
DEMAND - Computes Engine Demand 


AREAF - Computes Fy where Ay = FUWon 


Afterbody Module 

ABINPT - Reads Nozzle Input Deck 

DRAGAB - Monitors Afterbody Drag Caiculations 
DINTFR - Computes No@zle Interference Drag 
DBGATR - Monitors Boattail Dray Computation 
DBASER - Computes Base Drag 


DBTTB - Computes C, = £(M, PS) 
Der 


CDCGNV = Monitors C-D Nozzle Boattail Computation 


PLUGMX - Computes Boattail Drag for Plug Nozzle, 
Mixed Flow 


PLUGNM - Computes Boattail Drag for Plug Nozzle, 
Non-mixed Flow 


DRG - Computes Boattail Drag for C-D Nozzle 


SUBBT =- Computes Subsonic Cy for C-D Nozzle 
BT 


SUPET - Computes Supersonic Cy for C-D Nozzle 
BT 


CABTAB - Computes C as Function of A,,/A 
Dar 10° 9 


AEXHST = Computes Nozzle Exit Area 
Engine Module 
PERF <- Computes Gross Thrust 


CHKRP - Computes Thermodynamic Properties at 
Throat 
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GTHRST - Computes Performance of Propulsive 


Nozzles 
*) HOFT - Computes Enthalpy Given Temperature 
Ee TOFH -~- Computes Tempcrature Given Enthalpy 


PR@FH ~ Computes Relative Pressure Given Enthalpy 


- HGFPR - Computes Enthalpy Given Relative Pressurv 
C@FH - Computes Critical Velocity Given Total 
Enthalpy 
C@FHT - Computes Critical Velocity Given Static 
6 Enthalpy 


Miscellaneous Programs 
ST@RE - Stores results of a performance point 
€ @UTGS - Outputs an engine performance matrix 
WARNING - Prints warning message heading 
FIXDIM - Determine dimensions of A8 and A9 input 
« ATM@S - Computes atmospheric data 
DATA - Transfers and output input card images \ 
Tabular Input Programs 
€ TABLL - Inputs table format 1 tables 
TABL2 - Inputs table format 2Q tables 
TABL22 - Inputs table format 2K tables 
¢ TABL3 - Inputs table format 3 tables 
PLACIN - Reads a tabular array 


K@UNT - Counts non-blank input fields 


ro 


Linear Interpolation Programs 
TABUL - Interpolates F = £ (x) 


TABU2 - Interpolates F = f(x,y) 


«i 


vs 


Interpolates F = f (x,y,z) 
Finds Interpolation Interval 
Computes Interpolation Factor 


Perform Linear Interpolation Calculation 
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METH@D: 


USAGE: 


SUBPR@GRANS : 


SECTION I1 


DESCRIPTIONS OF SUBROUTINES AND LISTINGS 


PROGRAM TEM-333 


Main Program TEM-333 monitors the overall 
program logic. 


All supporting data is input, then the program 
reads performance points one by onc, computes 
and stores the results, and outputs an array 
of results whenever Mach or altitude changes. 


PROGRAM TEM333 (INPUT, @UTPUT, TAPL5, TAPE6 = 
4SUTPUT, TAPE] = INPUT) 


!NPUT or TAPE] System Input 

MUTPUT or TAPE6 System output 

TAPES Card Images of Input 
TAPE7 MARK II output 
ABINPT DATA FIXDIM 

PUTGP STORE TABIN 


TABL1 TABL2 
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: 


take te 


ON watt 5 RO 
(CALL:DATA] 


READ and WRITE 
TITLE 


READ ENGINE BLEED TABLE 
{CALL TABL 2} 


READ REFERENCE NOZZLE DRAG 
TABLE 


{CALL TABL1] 


READ INLET INPUT DATA DECK 


(CALL TABIN] 
a aH 
a J 
READ AFTERBODY INPUT DATA DECK 
{CALL ABINPT)} 
i 2 
| READ PUNCH CODE, NO. OF ENGINES, 
NO. CARDS/POINT, DATAS, AND DATA9 
> 
SET DIMENSION FLAGS FOR Ag AND 
FOR Ag INPUT FIELDS 
(CALL FIXOIM] & 


Figure 1: FLOW CHART FOR PROGRAM TEM 333 


PRINT OUPUT 
ARRAY FOR LAST 
* MACH — ALT 

{CALL OUTGO} 


NO 


ALT, FN, WET, W, Pp/P,. 
AB, A9, CFG, PS 


RINT OUTPUT 
ARRAY FOR LAST 


READ PERFORMANCE 
POINT DATA 
CARD IMAGE 


CONVERT A8 and AS 
INPUT FIELDS TO AS 

AND AQ 80 INCHES, IF 
REQUIRED 


Figure 1: FLOW CHART FOR PROGRAM TEM 333 (Cont) 


® 


YES 


NO ALT-MACH POINTS, 
NALT=1 


SAVE ALTITUDE 
[CALL OUTGO) 


La 


INCREMENT 
NO. OF ALT-MACH 
POINTS 

NALT NAL}+i 


ye 


THE INPUT POINT 
[CALL COMPUTE] 


COMPUTE INSTALLED ENGINE PERFORMANCE FOR 


No 


STORE RESULTS IN OUTPUT ARRAY FOR 
LATER OUTPUT. NALTTH OUTPUT COLUMN 


{CALL STORE] 
NEXT CASE 


(10 ) “f 


Figure 1; FLOW CHART FOR PROGRAM TEM 333 (Concluded) 
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USAGE: 


SUBPR@GRAMS : 


SUBR@UTINE COMPUTE 


SUBR@UTINE C@MPUTE controls the calculation 
of the installed performance for a single 
input case. 


Using the input performance data C@MPUTE first 
calculates atmospheric data and the reference 
recovery factor (MIL STD 5008 B). Gross thrust 
is computed for the input Wer Wor and Po/Po: 


The inlet module is called to find inlet drag 
and recovery factor. Then a second gross thrust 
calculation is made using values corrected by 

the recovery ratio. The input. array to the 
afterbody module is loaded and the module 

called to compute afterbody drags. Reference 
afterbody drag is obtained by interpolation. 

Now installed thrust is calculated by combining 
the appropriate terms and the total calculation 
is complete. 


CALL C@MPUTE 
All inputs are in labeled CQMM@N blocks. 
Input variables are: 

Mach, ALT, PS, P,/Po: Ag, Ag, Pye W a 
SFC, Wo , Table 2.1, W 
Cc 
and Table 2.2, C 


BLD = f (Mach, ALT), 
D 


= £(‘lach). 
AB por 
All of these input variables are direct input 
to the program or immediate consequences of 
the input, e.g. 


SFC = Wp/F yy: 


ATM@S DRAGAB INLDRAG 
PERF SQRT TABU] 


TABU2 


13 


~= 


) 
! 
! 
' 
l 


Figure 2: 
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Pomp {1 + .2 Mach?) 
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The Temp (1 +.2 Mach?) 
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FLOW CHART FOR SUBROUTINE COMPUTE 
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Mila, = 1. -.076 (Mach - 1.)1-38 


(CALL PERF) 


Camgute Om and 
RE « 1@acn Pro Ww) 
[CALL INLOAAG) 


c 
' 
' 
' 
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Social 


STORE FOR OUTPUT, ADDITIONAL 
WLET DATA 


De, “Deiy’ “Opyr’ “Owe 
* VA, AeE/ a. Ada. 


t Correct ‘ 
' fer ‘ 
’ Resovery Q 
em 


Figure 2: FLOW CHART FOR SUBROUTINE COMPUTE (Cont) 
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{ af 

: Compute 

’ Abateaty 
y Performance 
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Figure 2: 


Compute F, = Ls ene, Perens, Ty, Pama! 


(CALL PERF) 


COMPUTE AFTERBODY DRAGS 


“Der, Og ase, “Oinr, AB 
(CALL ORAGAS} 


FLOW CHART FOR SUBROUTINE COMPUTE (Concluded) 
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SUBJECT: 


PURPQSE: 


METH@D: 


USAGE: 


SUBROUTINES 
CALLED: 


SUBRQUTINE TABIN 


SUBR@UTIN. TABIN reads the Inlet performance 
data deck, Cards 3.1 to 3.2.12. 


The input deck is read. Table 2A is enriched 
as required to evaluate that table. If inlet 
capture area is not input it is computed. 
CALL TABIN 

All COMMON blocks in this routine are filled 
by calling this subprogram. 

SIZINL TABL1 TABL2 


TABL22 TABU] 
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USAGE: 


Output 
Output 


Input from 


SUBR@UTINE INLDRAG Controls the Calculation of 
Inlet Performance. 


After initialization M, is interpolated and 


compared with the minimum Mach number for 
which inlet drag will be computed and with 
the start Mach number to determine which inlet 
mode to use. Low speed mode is computed in 
INLDRAG: Subroutine AIRBYP computes external- 
compression mode performance, and subroutine 
AIRSPL computes mixed-compression mode. In 
the two latter cases inlet drag coefficients 
are interpolated and inlet drag determined. 
The outlet array is loaded and subroutine is 
complete. 


CALL INLDRAG (XM, P, WE, WS, DINLET, RI) 


XM - Free-stream Mach number 

P - Total free-stream pressure LB/FT* 

WE - Engine corrected airflow, primary 
LB/Sec 

ws - Engine corrected airflow, secondary 
LB/Sec 


DINLET - Inlet Drag LB 
RF - Inlet recovery factor 
COMMZN 


XMMIN - Minimum Mach for inlet drag caiculation 


AC ~- Capture area - £t? 
Inlet Tables 1, 2B, 3, 4 and 5 
Output in CZMMN 
SUTINL (10) 
SUBRGUTINES CALLED: AIRBYP AIRSPL AREAF 
TABU1 TABU2 
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es 


Mach, = f (Mach ) from Table 1; 


Apnea = AREAF (Mach,) 


Mixed - 


One Compression 
mode 


“Om ~ 20 
“gp, = 00 
Opyp = 0.0 
Ogio ~ 20 
Ao,/Ac = 0.0 
AolAc = 0.0 
Ao/Ac = 0.0 
RF = ¢ (Mach,) 
from Table 2B 


= f(Ao,/Ac, Mach,) from Table 3 


COgip =f (Aog/Ac, Mach,) from Table 4 


20 i 
Vv 
cp nD _~ 


“Ogyp ~ f1A0gyp Mach) from Table 5 
Ac 


Cone * COsp, *©Dgi5 *COpve_, 
0, “pF Mach? (1+ Tt Maoh?) ‘ee 


DINL = Qe Ac CDiyy 


OUTINLIT) = Cp... OUTINLTB) = Ao,7Ac 
OUTINL(2)= Cog, OUTINLIG) = Oi, 
OUTINLI3)=Cp,,,  OUTINL(7) = Ao/Ac 


OUTINL(4) = Ao,/Ac 


RETURN 


Figure 3; FLOW CHART FOR SUBROUTINE INLDRAG 
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SUBJECT: SUBR@UTINE SIZINL 
Pi'RPGSE: SUBRGUTINE SIZINL computes inlet capture area. 
USAGE: CALL SIZINL 


Inputs in C@MM@N/IDESIGN/ design Mach and 
airflow. 


Tables 1, 2B and 2C. 
Output: AC Et 


SUBJECT: FUNCTIQN AREAF 


Vad 
013858 (14 YS! mach2)? 7 1) 


PURP@SE: AREAF = : —~ 
Mach \/ Y 
USAGE: A = AREAF (XMACH) 
XMACH local Mach number 
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= “S' = 7 ae SS eT 
pees od 
~ Tonsil Ceenndines iieiieiane tetas ate ne * sage wena 


U SUBROUTINE ; 
8 siZine I 
Méchones = f (Mathes) from Table 1 
Proto = f (Machpes) from Table 2B 
Ao/Ac =f (Machpgc) from Table 2C ) 
Ao = AREAF (Macho, .<) We PJ, 
a mean aes a=» 
| FUNCTION 
AREAF | 


+1 
AREAF =  .013658(1+ = i Mach?) 2 '7-1) 


Mach VY 


RETURN 


Figure 4: FLOW CHART FOR SUBROUTINE SIZINL AND FUNCTION AREAF 
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SUBJECT: 


PURP@SE: 


METH@D: 


USAGE: 


SUBPR@GRAMS : 


SUBRGUTINE AIRBYP 


SUBRGUTINE AIRBYP Computes recovery and mass 
flow ratios for external compression mode inlets. 


AIRBYP matches the recovery and mass flow from 
Table 2A to that obtained from function DEMAND. 
Table 2A was enriched to 41 points per array 
when input by SUBR@UTINE TABIN. AIRBYP first 
interpolates in this table to develop a curve 
of RF =f (A/A,) at the input Mach number. 


The 4lst point of this table represents the 
choked flow point, in Maximum (Ay/AG) - DEMAND 


is called to compute (AQ /A,.) for the RF at the 
choked point. This A/AG is compared to the ASA 
of the choked point. If the Ad/Ac from DEMAND 

is greater than the maximum Ad/Ag then the inlet 


is choked and RF must be found such that 

DEMAND (RF) = Ad/Ac . A Newton-Raphson iteration 
MAX 

is used to solve this case. Otherwise the iniet 

is not choked and a solution for A,/A, exists 

less than the maximum. A Regula Fist iteration 

is used to match the recovery and mass flow. 


In either case buzz and distortion limits are 
checked and the other mass flow ratios computed. 


CALL AIRBYP : ' 
Inputs in COMMON include Mach number, air flow, 
capture area and inlet tables. Other COMMON 
blocks are output. 
DEMAND FACINT IFIND 


TABU2 


TABU1 WARNING 


START 
Y) 
GUILD RF = f (Adp/Ac) CURVE FOR INPUT MACH NUMBER. LINEAR 
INTERPOLATION ON MACH IS APPLIED POINT BY POINT TO THE 
ENRICHED TABLE 2A. THE 41ST POINT REPRESENTS THE 
CHOKED CONDITION. > 
No Old Curve 
Is 
Still Good & 
Yo 3 
J 


INTERPOLATE FOR MACH IN ENRICHED TABLE 2A THUSLY: 
FOR K=#= 1,41 


INTERPOLATE FOR MACH IN ENRICHED TABLE 2A THUSLY : 
FOR K#=1,41 P(K) = XMFAC (PTAB (K, | + 1) — PTAB (K, 1))+PTAB (K, 1) 
A(K) = XMFAC (ATAB (K, 1+ 1) — ATAB (K, 1)) + ATA (K, 1) 
WHERE 11S LGWER SOUND OF MACH INTERNAL AND XMFAC = Mach — Mach ;), 


Machy) + 1) — Mach (1) 
[IFIND and FACINT] 


CHECK FOR CHOKED FLOW CONDITION 


4 FUNCTION DEMAND COMPUTES]| 
i Ao/Ac = (RF) 


FC = C-A2 =P on Gb OD aw © oe ow ow a ow a ee 


Az = DEMAND (P(41)) 


C<Ay “CHOKED 
FLOW” 


CPA, 


Figure 5: FLOW CHART FOR SUBROUTINE AIRBYP 


REGULA FALS/ 
ITERATION 


TO FIND 
INTERSECTION 


NON-CHOKED FLOW 


R1 = f(Ao/Ac) from Curve 
A2 = DEMAND (RF) 
F = Ao/Ac — A2 


F-FO 


Figure 5. FLOW CHARTS FOR SUBROUTINE AIRBYP (Cont) 


NEWTON-RAPHSON 4 
ITERATION TO FIND 
INTERSECTION 
CHOKED FLOW ; 
3 
je) Ao/Ac = A(41) 
F = FC 
| SECOND STARTING POINT ” P 
RF = P(41) { 
Xq = RF/2 
Fa = Ao/Ac — DEMAND (RF) 
¥ ad I = 0 ; 
X,, = RF - F (RF - Xo) 
‘wis ou F -. 7 2 } 


Fo =F 
RF = X, 
F = Ao/Ac — DEMAND (RF) 


> B 


Figure 5: FLOW CHART FOR SUBROUTINE AIRBYP (Cont) 


“ITERATION 7 Last Values 
| RF Used | 
| eee Tas To a sa awd 
Check for 
D - Buzzand 
Distortion 
Limits 


Mach < TABs px (1) 


| Minimum Mach | Mach #TABony (1) 


in Table 20 
isMinimum for | Buzz = f (Mach), fr..n 
L_ a= Ga» en oe J 


Ao/Ac < BUZZ 


“Ao/Ac BELOW 
BUZZ LIMIT” 


” Ao/Ac : BUZZ 


Ao/Ac # BUZZ 


DISTOR = ¢ (Mach,) from 


Ao/Ac >DISTOR “Ao/Ac Exceeds 


= = f(Aog/Ac, Mach,) from Table 7 


= #(Ao/Ac, Mech.) from Table 6A mete 


Figure §: FLOW CHART FOR SUBROUTINE AIRBYP (Concluded) 
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SUBJECT: 
PURPQSE: 


METH@D: 


USAGE: 


SUBPR@GRAMS : 


SUBRGUTINE AIRSPL 


SUBR@UTINE AIRSPL computes recovery and mass 
flow ratios for mixed-compression mode inlets. 


The various mass flow ratios are computed from 
the input airflow and by interpolation in 
Tables 6B and 2C. Recovery comes from Table 2B. 
If bypass mass flow is negative; recovery is 
recomputed so that the bypass mass flow is zero. 
A message stating that the inlet is undersized 
results. 


CALL AIRSPi’ 
Inputs in C@MMPN include Mach number, air flow, 


capture area, and inlet tables. Other CQ@MM@N 
blocks are output. 


TABU1] WARNING 


5 


= f (Mach) from Table 68 
a 


Ms 
Ac 
Ao 
kK 


= {(Mach,) from Table 2C 


AOL = Ac/dc + Aog/Ac 
Ac 
RF = ¢ (Mech,) from Table 28 


Ade Area W, RF 
Ac Ac 


ad 


— = Ao/Ac - Adg/ Ac 


Ao 
BY 
Ac 7 90 


RF = Ao/Ac = i 
Aree W.. 


RETURN 


Figure 6: FLOW CHART FOR SUBROUTINE AIRSPL 
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SUBJECT: SUBRGUTINE DRG 


PURPQSE: SUBR@UTINE DRG computes bottail drag for 
convergent-divergent nozzle. 


METH@D: The built in nozzle data tables are defined in 

: the DATA statements of this routine. Boattail 
angle is computed; then Mach is interrogated. 
If Mach is less than .9 the subsonic procedure, 
SUBBT, is used. Similarly if Mach is greater 

7 than or equal to 1.0 the supersonic procedure ' 

is used. If Mach is in the range from .9 

to 1.0, then the subsonic procedure is used 

with Mach equal to .9 and the supersonic 

procedure with Mach equal to 1.0. Finai results 

: are obtained by interpolating between the two 

. sets of results using Mach as the independent 
variable. This provides a smooth transition 

between the two Mach regimes. 


| USAGE: CALL SUBR@UTINE (DMDG, DBDJ, LDJ, EMF, EMJD, 
PIPF, THETA, GAM, ICQN) 
: | DMDJ Ratio of maximum to jet diameter 
| DBDJ Ratio of base to jet diameter 
‘ LDJ . Ratio of length to jet diameter 
| EMF Freestream Mach number 
EMJD Jet design Mach number 
PTPF Ratio of Jet total to free-stream x 
A static pressure 
| THETA Boattail trailing edge angle - degrees 
| GAM Ratio of specific heats 
ICON = 0 circular are boattail 


1 conical boattail 


Output in CYMM@N/A/ 


SUBPR@GRAMS : ATAN SUBBT SUPBT 


oan aM 


START 


free weaeweaw =n oe On an oe a 


i Define Built-in Dets 5 
fa cs cen ee tare ‘an om om oll 


Toh = (2a _ Pas. 

Dg 09 
2 (L/D9) 

Be = Tan7! (Tan (61) 


Mach 1.0 . 
Using Supersonic Method 
[CALL SUPBT) 


Compute “Dgy Og... and°0,,, 
Using Subsonic Method 
(CALL SUBBT) 


Compute “Der, “Dae, and “Dos, Using the Subsonic 
method with Mach Number Equal to .9 
(CALL SUBBT] 


e- 
= 


Compute °Dgy. “OB ese. and “Dror, Using the Supersonic 


Method with Mach Number Equa! to 1.0 


(CALL SUPB8T} 


-.9) (Op - Dp, ) 


—Dpase * “Opens, + Mach = 01 (Op eee, — Opage,) 
1.0-.9 


“Dro 7 “Orot, * {Mech - 9) (CD yo, — Or o,,) 
x 


1,0 — 
Figure 7; FLOW CHART FOR SUBROUTINE DRG 
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SUBJECT: SUBRQUTINE CDCQNV 


| 
PURP@SE: SUBRGUTINE CDCZNV moniters the convergent- 
divergent nozzle boattail drag computation. ) 
METH@D: CDC@NV computes the required input values for 


SUBR@UTINE DRG from the TEM~333 inputs. If 

the supersonic procedure is to be used }.e., 
Mach 21.0, these inputs include the design 

Mach number of the nozzle which is solved F) 
for by iteration. DRG is called; then boattail 

drag computed from the CD obtained. 


Note: The set DRG, SUBBT, and SUPBT were in 
existence prior to TEM-333, hence the incorpor- ; 
ation of this routine into the logic. 3. 


USAGE: CALL CDCQYNV 


Input through C@MMZN/ABIN/ and C@MMZN/AMAX/ 

Output includes D,. for interference routine 

Flags PLUG = 0.0; BITAB = 0.0 causes this a 
subroutine to be called. 


SUBPR@GRAMS : AEXHST DRG SQRT WARNING 


Diemeter for Interference Drag 
0, = OF 


7.9 


Solution for Nozzle Design Mach Number 
Newton—Rapheon Technique 


Find Pa/P, Ratio Such that Ag/Ag Equals 
that of Input Performance Point 


(FUNCTION AEXHST) 


1 = 0 €= 0008 
X1 = Pa/Po 


XQ = (1.1) Py/P, 


Figure 8: FLOW CHART FOR SUBROUTINE COCONV 


to Insure ' 
Convergence 4 
of Low ‘ 


PaO, 


P ) 
Fy = AEXHST (X4) — Ag/Ag 
2 
Fars AEXHST (Xo) - Ag/Ag 
F, =F 
1 2 
a 
F, #F, 
1<5 

| 1:5 Pig = 1.09 ~ (.02)! 

4 


Pa/Pon = Max Me/Pop, Pix! | 


F = AEHST "7/Pop)—Ag/Ag | 
| 


a 


< © 


Figure 8: © FLOW CHART FOR SUBROUTINE COCONV (Cont) 
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twee 


= f (Dy oa L_ Mech, Mp Pg /Po, 6, y, ICON) 
B. 65. Oy 


(CALL ORG) 


RETURN 


Figure 8: | FLOW CHART FOR SUBROUTINE COCONV (Concluded) 
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SUBJECT: SUBRGUTINE PLUGMX 
PURPQ@SE: . SUBRGUTINE PLUGMX computes boattail drag for a 

plug nozzle, mixed flow. ) 
METHQD : Boattail angle is computed from the maximum 


diameter, exhaust diameter, and plug length. 

Then Mach is checked and either the subsunic 

or supersonic method applied. Input tables 

4.2.2 and 4.2.3 support the supersonic J 
procedure. Boattail drag is returned. 


USAGE: CALL SUBR@UTINE PLUGMX 
Input from C@MM@N/ABIN/ and C@MM@N/AMAX/ 
Output includes De for interference routine 3 
Flags set such that: PLUG = 1.0 
FLWMIX = 0.0 
BTITAB = 0.0 : 


Will cause this routine to be called 
SUBPRZGRAMS : AEXHST ATAN SQRT 


TABU] TABU2 


4 ete mett etree meee oe Sree om eer nt res te AAA CLE CNS CELE RAR teh ie ko asim oo meena eee AONE he ROO INEM ed Pe 8 eaten 


4Acx 


> 0.0 


<1.0 71.0 


f (Mach, 8) trom Table 4.2.2 


Cas 1.4 Tan (0) (1- Age) 
B Amex 


f (Pg/Po) from Table 4.2.3 
= 10- (1--Cvy) 
P HEP 
c 
Ce = %v, - O Amax 
9; 


(CE, — Ceo) Fo, 


Mech 1-53 


Pat * Cog O- Amex 


RETURN 


Figure 9: FLOW CHART FOR SUBROUTINE PLUGMX 
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SUBJECT: SUBR@UTINE PLUGNM 
PURPQSE: SUBRQ“UTINE PLUGNM computes boattail drag for 

a plug nozzle, non-mixed flow. a 
METH@D: Boattail angle is computed from the maximum 


diameter, the exhaust diameter and the plug 

length. Then Mach is checked and eithcr the 

subsonic or supersonic method applied. Input 

tables 4.2.2 and 4.2.3 support the supersonic 2 
procedure. 


Boattail drag is returned. 


USAGE: CALL SUBRGUTINE PLUGNM P 
Input from C@MM@N/ABIN/ and C@MM@N/AMAX) 
Output includes De for interference routine 
: Flags set such that: PLUG = 1.0 
FLWMIX = 0.0 ; 
BTTAB = 0.0 
Will cause this routine to be called 
SUBPR@GRAMS : AEXHST ATAN SQRT : 
TABUL TABU2 
‘ 


a 
= Dy 


f (Mach, 8) frorn Table 4.2.2 
f (Pg/Po) from Tabie 4.2.3 


f (Pyg/Po) from Table 4.2.3 


= Cc 
0.5 ( Vp + Cy) 


10 - (1- Cy.) 
y Buep* ie 
(CEg — CE p90! (Faint Fine.) 


Figure 10: FLOW CHART FOR SUBROUTINE PLUGNM 
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SUBJECT: SUBR@UTINE ABINPT 


PURP@SE: SUBRGZUTINE ABINPT reads the nozzle performance 
data deck, cards 4.1 to 4.2.6. 


METH@D: The input deck is read. The thickness of the 
annular base is doubled to account for top and 
bottom when compared to the total jet diameter. 
If Alo was input the table of Mach, vs. Ag/Ag 


needed by subroutine CABTAB and the P 9/P amp 


lower thrust are computed. Tables not indicated 
by the input flags are not read. 


USAGE: CALL ABINPT 


All C@MMZN blocks present are defined by 
Calling this subprogram. 


SUBPR@GRAMS : TABL1 TABL2 TABL3 
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SUBJECT: SUBRQUTINE SUBBT 
PURPQSE: SUBRGUTINE SUBBT computes boattail drag for 
a convergent-divergent nozzle when Mach < .9. 5 
: , 2 ae 
METHQD: The quantity D, Ps - Dy) is tested. If it is 


greater than 0.25 boattail drag is determined 
directly from built in tables XI and XII. If 
not, an emperical method to develop base and 
boattail drag is used. 


USAGE: CALL SUBBT (DMNDJ, DBDJ, EMF, PTPF, THETA, 
IC@N, GAM) 
DMDJ Ratio of maximum to jet diameter ' 
DBDJ Ratio of base to jet diameter 
EMF Free-stream Mach number 
PTPF Ratio of jet total to free-stream : 
static pressure 
THETA Boattail trailing edge angle - degrees 
GAM Ratio of specific heats 
ICGN | = Q circular arc boattail : 
= 1 conical boattail 


SUBPR@GRAMS : CABTAB TABU] TABU2 


Cog, = £18, (0B/D))/(OnvD)) from Table IV 
ACy 


= (8) from Table V 


Pt Pf: .001 


ACh/C,= t (8) trom’Table VIII 


ACA/C.=f (8, Pe/Pe) from 
op Taek vi 
AU A Cp,)=#(Py/Pe) trom Table tit 


Figure 11: FLOW CHART FOR SUBROUTINE SUBBT 
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B Cp =F (PPoo) 
table X11-A 


from table XIIB 


Pa/P_ = f (P MPooltrom table Xill 
= @2 (Pa/P.. =1,0) “1, 
“on, a; (ay 1 ) . 
RANG ees ames 
ge 


Figure 11: FLOW CHART FOR SUBROUTINE SUBBT (Cont) 
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SUBJECT: 


PURPQSE: 


METHQD: 


USAGE: 


SUBPRYGRAM: 


SUBRQGUTINE SUPBT 


SUBRGUTINE SUPBT computes boattail and base drag 
for a convergent-divergent nozzle when Mach > 1.0. 


Boattail drag is computed using Table I, then 
if a base has been defined base drag is computed 


and added to compute a total drag. 


CALL SUPBT (DMDJ, DBDJ, EMF, EMJD, PTPF, GAM, 


TANB) 


DMDJ 


ALQG 
SQRT 


Ratio of maximum of jet diameter 
Ratio of base to jet diameter 
Free-stream Mach number 

Jet design Mach number 


Ratio of jet total to free-stream 
static pressure 


Ratio of specific heats 

Tangent of boattail chord angle 
CABTAB EXP 
TABU1 


be K+] 


“Ogr = (Ag/A ro) 
(CALL CABTAB) 


PP; =(1+ p= Mogg’)! 77" 
205 .. 1M 
Cpa pes 


C= (.37 Mach + .62) =Des 
a 
( 0; ) PP; 


B=f(C) from table !I 


K = EMF 
Pa/Pe = £(K) “rom table II! 
Cy. = 72 (Pye - 1.) 
yMech* /0 
(oa) 


2 


Pa’ = f (K) from table ti! 


Pa/Pe = Pa/Pe B 
Meche 815 - 1.15 ALOGIK) 


RETURN 


Cp +C 
0” “Dey Ogs 


Figure 12: FLOW CHART FOR SUBROUTINE SUPBT 
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SUBJECT: 
PURPQSE: 
| METHQD: 
| USAGE: 
| 
ABIN 
1 AES 
2 As 
3 AELB 
4 alg 
5 Pa/Py 
6 Py 9/Po 
7 FGIP 
8 FGIF 
9 PC 
10  GAMMAP 


ene rn 


SUBR@UTINE DRAGAB 


SUBRGUTINE DRABAG monitors the afterbody drag 
calculation. 


The input array is transferred to the afterbody 
CQMMZ<N block ABINP. A check is made on Mach 
number. If Mach is equal to zero drag is set 
equal to zero and no further afterbody calculations 
performed. If Mach is greater than zero, boattail 
drag, nozzle interference, and base drag routines 
are called to compute the various drags. The 

total airplane drag is returned for nozzle 
interference and base drag so these are converted 
to one engine values. The output array is 

loaded with the various C,'s based upon area 

AMAX (drag may have been “computed using a 
different area). 


CALL DRAGAB (ABIN, ABQUT) 

ABIN (40) - Input array 

ABZUT (10) - Output array 

The input array is defined as follows: 


Primary nozzle effective area - in? 


Primary nozzle geometric thrust area - in? 


Fan nozzle effective area - in? 
Fan nozzle geometric thrust area - in? 


Primary nozzle pressure ratio 
Fan nozzle pressure ratio 


Primary nozzla ideal thrust (actual mass flow) lb 
Fan nozzle ideal thrust (actual mass flow) lb 
Power setting 


Primary flow specific heat ratio 


ABIN (continued) 


11  GAMMAF 
12 XM 

13 PAMB 
14 Q 

15 Ag 

16 §=Al19 

17. CFG 

18  DHEP 
19 DHEF 
20 PLUG 
21 FLWMIX 
22 DMAX 
23. THICKB 
24 ZLPLUG 
25  ZULBTAL 
26)6=—Cl S 

27 ~=ABASE 
28 BTTAB 


32-40 
SUBPR@GRAMS : 


Fan flow specific heat ratio 
Free-stream Mach number 

Ambient static pressure psia 
Free-stream dynamic pressure psia 
Primary nozzle exit area in? 

2 


Fan nozzle exit area in 


Nozzle internal performance coefficient 


Nozzle type designation 

Mixed flow flag 

Maximum diameter, in. 

Total thickness of annular base, in. 
Overall plug length, in. 

Nozzle boattail length, in. 

Nozzle centerline spacing, in. 

Base flow area, in’. 


Boattail table flag 


Number of nozzle interference spacings 


Not defined 


DBASER DB@ATR DINTFR 


LOAD LOCAL VARIABLE 
ARRAY FROM INPUT 
ARRAY 


ee No Is Yes ser ate se 
AX*"4 Onsax Mech OUTPUTSTO | 
= 0? ZERO _J 
Compute boattail 
Drag, Dgy and 8 
[CALL DBOATR } Cag = 0.0 
> 00 
“Der 
= 0.0 
Compute total A/P “Dwi 


Nozzle interference Drag, 
ON! A/P 
(CALL DINTFR] 


= Q0 
“Dpase 


B= 0-0 


Compute total A/P 
Beee Drag, D 
rag Base a sp 


[CALL Dgase a! 


Ont ~~ Nt wp (jee) 


LOAD OUTPUT ARRAY 
Dag ~ Spy * On, * Opase 


= Dgq/(Q: AMAX; 
Cy, Ont" AMAX) 
Coe age * Oase/(Q: AMAX) 


Figure 13: FLOW CHART FOR SUBROUTINE DRAGAB 
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SUBJECT: 
PURP@SE: 


METH@D: 


USAGE: 


SUBPR@GRAMS : 


SUBRGUTINE DBYATR 


SUBR@UTINE DBZATR monitors the calculation of 
boattail drag. 


The input flags PLUG, FLWMIX, and BTTAB are 
interrogated to determine which boattail ; 
method is desired and the appropriate routine °. 
selected. 


CAL DBSATR 
Input through C@MM@N/ABIN/ 
CDCZNV DBTTB PLUGMIX PLUGNM 


Ne 


5 
a 


 < ~~ Tie St, og FOP Pa Ae Ome * 


we 
™ a 


é 
START 
N <0 
Oar =0.0 
q 
uel OB is interpolated 
option? from table 4.2.1 
(CALL OBTTB) 
y No 
i“ 
No CALL CONVERGENT 
NOZZLE MODULE 
(CALL COCONV]) 
ad 
C Yes 
Yes CALL MIXED FLOW 
PLUG NOZZLE 
( MODULE 
(CALL PLUGMX] 
No 


CALL NON-MIXED 
FLOW PLUG 

NOZZLE MODULE 
(CALL PLUGNM) 


am 


RETURN 


Figure 14: FLOW CHART FOR SUBROUTINE DBOATR 
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SUBJECT: 


PURP@SE: 


METH@D: 


USAGE: 


SUBPR@GRAMS : 


SUBRGZUTINE CABTAB 


SUBR@UTINE CABTAB computes afterbody drag as a 
function of Ag/Ayg: Po /P amp’ and Mach. 


The drag coefficient is input in tabular form 
in input Table 4.2.4. Mach and A,/A39 are known; 


the suprogram computes Po /P MB after interpolating 
Mach, = £(Ag/Ag) from a table generated by the 


program. The table senerated in the input 
routine ABINP is developed as follows: 


For Mach, = 2.8, 2.6, 2.4 ..., 1.0 compute 
y+1 
sa race -y+] 
2(y-1) Crs 
= yt z, y-1 2\ 2 (y~1) 
Ag/a, = (S+) Machy (1+%5+ Mach,”) 


The resulting array of Ag/Ag is strictly 


increasing; hence serves as an independent 
variable array with Mach, as the dependent 


9 
array. The subroutine also checks P,/P against 
1 lower bound value of PY/Po = 1/( 


LB 
The C obtained from Table 4.2.4 is based on 


Das 
an area Axo whereas the program uses Auax (most 


of the time these two inputs will have the same 
value), consequently the Cy is ratioed by Ay o/Awax 


before returning to the calling program. 
CALL CABTAB (CAB, AM) 

Output CAB - Afterbody drag coefficient 
AM Free:‘stream Mach number 


Input from C@MM@N/ABINP/,/AMAX/,/AREA10/,/P8PAL/, 
/XMACH9/, and /ABTAB6/ 


TABU] TABU3 
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2 
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the 


ary 


Ne 


a ees enti gene 


fN 


bf 


<(Pg/Pp) LB 
Po/P a : «1 


> (Pos) Lp 


Ye <1.0 


Mach = tiAg/Ag) 


from internelly 


! Note 1 (Palos = —_,- 


x | — 


a H (a) ot 
Wego ee ee ee 


Internal table, For Macha = 2.8, 2.6, 2.4, ..., 1.0 
r tL 


1 helMg -(221) 2(7-1) 


i MglAg array will be strictly increasing 
Figure 15: FLOW CHART FOR SUBROUTINE CABTAB 
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SUBJECT: SUBROUTINE DINTFR 
: PURP@SE: SUBROUTINE DINTFR computes nozzle interference 
drag. 
METH@D: The ratio 8/D is interrogated to determine if 


it is greater than 1.0. If so, Cy is 


interpolated from Table 4.2.6 and total inter- 
\ ference drag computed. 


USAGE: CAL DINTFR 


Input through COM@WIN/ABIN/ except D which is 
input in CQ#MOQIN/DSPACG1. D is an output from 
' the boattail routines except when Cy) is 
BT 


interpolated from Table 4.2.1. In this latter 
case D = 1.0, then nozzle interference drag 
can be input only as a function of Mach and 
some constant S/D ratio where the ratio is 

| input instead of the nozzle centerline spacing 
for variable S. 


SUBPR®GRAMS : TABU2 


Note: Ocomputed by bosttall drag routine 


Figure 16: FLOW CHART FOR SUBROUTINE DINTFR 
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SUBJECT: 


PURPQSE: 


METH@D: 


USAGE: 


SUBPR@GRAMS : 


SUBR#UTINE DBASER 


SUBR#UTINE DBASER computes total airplane base 
drag. 


If a base area is input the base pressure 
ratio is interpolated from Table 4.2.5 and 
drag computed. 
CALL DBASER 

Input form C@MMPN/ABIN/ 


TABU1 


101 


es 


Pa = # (Mach) 


Pamb 


frorn table 4.2.5 


ve 


Pp 
B 

Opase ae = ) Pamb ABASE 
amb 


RETURN 


Note: Ag ASE 8 input in field 8 of card 4.1.1 


Figure 17: FLOW CHART FOR SUBROUTINE DBASER 
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_ SUBJECT: 


PURP@SE: 


METHQD : 


USAGE: 


rae SORE stance eee ane — 
= oo = s 
a 


SUBRGUTINE PIXDIM 


SUBRMUTINE FIXDIM determines the dimensions 
of the nozzle throat and nozzle exit inputs. 


The input flag is interrogated and three flags 
are returned indicating input type and 
dimension. 

CALL FIXDIM (FLAG, J, K, L) 


FLAG ~- Input code value 


Cy 
v 


0 area input, = 1 diameter input 


1 feet 


x 
i} 


0 inches, 


L = 0 ignore, 1 (Ag/Ag) input. 
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SUBJECT: SUBRGUTINE ATMZS 
PURPOSE: SUBR@UTINE ATMYS computes ambient pressure and 


temperature using empirical equations for the 
1962 U. S. Standard Atmosphere. 


METH@D: Altitude is tested to determine which region 
it is in, then the appropriate eqiations are 
applied. If altitude exceeds 100,000, an 
error is assumed and altitude will be divided 
by ten until it is within the required range. 

USAGE: CALL ATM@S (ALT, TAMB, PAMB) 

ALT ~ pressure altitude 
TAMB - ambient temperature 
PAMB - ambient pressure 


SUBPR@GRAMS : EXP WARNING 


110 


Tomb = 518.68 — 00956618 ALT 
ee 


Tomb 


Pamb = 21 10.22( 518.60 


T, = 380.97 
emb 36200 - ALT 


Pont = (2116.22) (.222408) ¢ 20010.7 


Tomb * 364.18 + 0008439 ALT 


30.4757 — .086643 T 
Pomb = 2116.22 « ame 


ALT EXCEEDS 


i204 ALT= .1ALT 
» ERROR 


| 
RETURN 


Figure 18: FLOW CHART FOR SUBROUTINE ATMOS 
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SUBJECT: 
PURP@SE: 


METH@D: 


USAGE: 


SUBPR@GRAMS : 


To re eee Neneeeeme ee 8 


SUBRGUTINE PERF 


SUBRGUTINE PERF controls the calculation of 


gross thrust. 


Mass flow, enthalpy, pressure, and relative 
pressure are developed for the nozzle throat. 
The thermodynamic properties at th2 throat are 
computed, then gross thrust computed from 
propulsive nozzle performance analysis. 


CALL PERF (WET, WZ, PT8@PA, TT, PAMB, FG, A8) 


Input WEFT - 
WE - 
PT8MPA- 
TT@ - 
PAMB - 
Output FG - 
A8 = 


CHKRP 


Fuel flow lb/hr 

Air Flow lb/sec 

Nozzle throat pressure ratio 
Total temperature degrees R 
Ambient pressure PSF 

Gross thrust lbs 


Nozzle throat (given) ft? 


GTHRST PROFH 


START 


FOAg = Wy, 
FOA, =0°0 


Ha = HOFT (T, , FOA2) 
Wg = Wo + Wy 
Hy = Hp + 18400. Ws Tg 
We 
Pg = (Pa/Pa) Pamb 


Prag ~ PROFH (Hg, FOAg) 


COMPUTE CHOKED THROAT 
PROPERTIES 


[CALL CHKRP} 


COMPUTE PROPULSIVE 
NOZZLE PERFORMANCE 


[CALL GTHRST) 


RETURN 


Figure 19: FLOW CHART FOR SUBROUTINE PERF 
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SUBJECT: SUBRGUTINE QSUTGS 


PURP@SE: SUBRGUTINE ZUTG®S prints a matrix of performance 
cutput. 5 
METH@D: The output for NALT power settings at a Mach- 


altitude combination has been stored in the 

array @UTP. A FORMAT is built to print the 

variable name on both the left and right of the 

output columns. The matrix is printed and any j 
indicated punching of MARK II data is done. 


USAGE: CALL QUTG# (NALT) 


NALT - Number of columns in output 
matrix = 19. ’ 


Input in CYMMON 


@UTP (45,10) output matrix 

PC@DE punched output flag 

DA, TE date, alphanumeric 

TITLE (8) run title 

MACH Mach number J 
Output | 

IWARN Error message indicator 

introduced 


SUBPR@GRAMS : MIN@ M@D 
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SUBJECT: 
PURP@SE: 


METH@D: 


USAGE: 


- — nen coprengern, 
oh LR OTC tp tet emer 
oH ete gern eran eney oe = Ft nee ener ee em ere ethene ot $ 
a 


SUBR@UTINE WARNING 


which the warning which follows applies. 
If no warning massage has been written since 
the last printing of an output matrix, a heading 
is printed. In any event, the case number, 
Mach, altitude, and power setting is printed. 
The calling subprogram will then print the 
appropriate message. 
CALL WARNING 
Input in CQ@MMON 

IWARN Error message flag 

CASE Case number 

ALT Altitude 

MACH Mach number 

PS Power Setting 


1 

SUBR@UTINE WARNING identifies the point to 3 
i 

| 

| 

| 

Output: STATUS Error flag for case. 
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SUBJECT: 


PURP@SE: 


METH@D : 


USAGE: 


SUBR@UTINE STPRE 


SUBRGUTINE ST@RE moves the results for a single 
case to the output matrix. 


The C@MMON block @UTS contains the output for 
the case just computed with the individual 
variables in exactly the order to be printed. 
This data is transferred to COMMON block ZUTP. 
CALL ST@RE (NALT) 


NALT - number of cases at current 
Mach altitude = 10 


@UTS (45) - output for naLT?# case. 


NVAR - number of ZUTS locations 
being used. 


Output OUTP (RT,10) Output 'matrix. 
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SUBJECT: 
PURP@SE: 


METH@D: 


USAGE: 


| 
| 
| 


SUBPR@GRAMS : 


SUBRGUTINE CHKRP 


SUBR@UTINE CHKRP calculates the thermodynamic 
properties at the throat and the total/static 
pressure ratio at the sonic condition. 


FUNCTIQN C@FH is called to compute sonic 
velocity. Then enthalpy is computed using 
the entry total enthalpy and the source 
velocity. Critical static temperature and 
relative pressure are computed as functions 
of enthalpy and fuel-air ratio. Pressure is 
derived from entry pressure and the relative 
pressure. 


CALL CHKRP (Pl, Hl, PR1, F@Al, PAM, TC, C, PC, 
P10PAM, PIPCC) 


Pl Entry total pressurg PSF 
Hl Entry total enthalpy, BTU/1b,, 
PR1 Entry relative pressure 
FQA1 Entry fuel-air ratio 
PAM Ambient static pressure, PSF 
TC Criticai static temperature, 
degrees R 
Cc Critical sonic velocity ft/sec 
PC Critical static pressure, PSF 
P10PAM P1/PAM 
PIPCC P1/PC 
C@FH PREFH T@FH 
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SUBJECT: 


PURPQSE: 


METHQD : 


USAGE: 


Output 


SUBRGUTINE GTHRST 


SUBRGUTINE GTHRST computes the performance of 
propulsive nozzles. 


Outlet conditions are computed, then nozzle 
ideal exit velocity. Nozzle throat static 
pressure is compared to ambient pressure 

to determine which temperature, pressure, 

and velocity to use in the calculation of 
nozzle throat area. Gross thrust is computed 
by the appropriate equation. 


CALL GTHRST (XN@ZZ, PR1, Pl, Hl, Wl, FQ@A1, 
PAM, CV, CD, TCS, PCS, VC, H2, T2S, V2, A, FG) | 


XN@ZZ - 1-0 for convergent-divergent nozzle; 
2.0 for convergent nozzle. 


P121 Nozzle inlet relative pressure 

Pl Nozzle inlet total pressure PSF 

Hl Nozzle inlet total enthalpy BTU/LB,, 

W1 Nozzle inlet mass flow LB/SEC 

F@A1 Nozzle inlet fuel-air ratio 

PAM Ambient static pressure PSF 

CV Nozzle velocity coefficient 

cD Discharge coefficient 

TCS Nozzle throat static temperature 
degrees R 

PCS Nozzle throat static pressure PSF 

vc Nozzle throat velocity ft/sec 

H2 Nozzle outlet static enthalpy BTH/LB,, 

T2S Nozzle outlet static temperature 


degrees R 


a, ata aan BEA 
ne + RS : 7a 


+ ee cement eiinwnee. ne 


; 
: Output (continued) i 
V2 Nozzle ideal exit velocity ft/sec 
A Nozzle throat area £t? : | 
FG Gross thrust lb 
SUBPRZGRAMS : AREA H@FPR SQRT 
TOFH 
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: 
por] 


PR aa (PR, iP ‘amb! 
+ 


H. = HOFPR (PR5, FOA,) 


Tog = TOFH (Hp, FOA,) 


4H= H,- He 


4H = 0°0 < a0 


A =0'0 4H 


IRROR = 1 


2 A+ emer nee TS te etn ES oN ese 


> 0.0 


cu 


Peg >Pamb 


A= AREA (W1, Tog, Pes, Vc, CD) 


Fee —- ~~ 


a er a 


RETURN 


Figure 20: FLOW CHART FOR SUBROUTINE GTHRST 
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SUBJECT: 


PURP@SE: 


METH@D: 


USAGE: 


t* 


os 


SUBRGUTINE DATA 


SUBRGUTINE DATA provides a gard image output 
of all data cards input to TEM-333. 


Data cards are read from TAPE] (INPUT), 
written on TAPE5 and OUTPUT. A card count 
is provided. Upon sensing an end-of-file 

of TAPE], TAPE5 is end-filed and rewound. 
All subsequent input is from TAPE5. The 
INPUT file is equivalent to TAPE1 to provide 
a logical number for the EGF test. 


CALL DATA 


INPUT file must be in desired position when 
subroutine is called. 
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SUBJECT: 


PURP@SE: 


METH@D: 


USAGE: 


SUBPRYGRAM: 


SUBRQ<UTINE PLACIN 


SUBRGUTINE PLACIN reads a tabular array. A 
count of points is provided for independent 
variables. 


If the number of fields input, NF, is equal to 
zero, then data is read until the maximum 
allowable number of cards has been read or 
until a 1 is input in column 72. A field by 
field count is made on the last card image so 
that thea total number of non-blank fields is 
returned in NF. If NF is positive, then NF 
values are read. All arrays are formatted 

ten field-seven digit. 


CALLPLACIN (A(1),NF,NCARD) 


A - array to be filled 


NF - number of fields to read or if zero 
on input, the number of fields read. 


NCARD - the maximum number of cards to be 
read if array is full. 
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SUBJECT: 


PURPQ@SE: 


METH@D: 


USAGE: 


ST@RAGE: 


Address 


Content 


SUBPR@GRAMS : 


ee niente A eS en A ta ee 


SUBRGUTINE TABL1 


SUBRGUTINE TABL1 inputs tables of the form 
F = £(x), table format 1. 


SUBR@UTINE PLACIN is called to read and count 
the independent array, then recalled to read 
the dependent array. 


CALL TABL1(T(1), MAX, N) 


= - Table to be read 
MAX - Maximum number of points in X irray. 
Output - N - Actual number of input points 


in X array. 


The table T is stored as follows: 


T (Max) T mMaxtl) 7(2 Max) 


N MAX 


T can be specified in the calling program by 
either: 


DIMENSION T (2MAX) or C@MMON/T/ TX (MAX) , TY (MAX) 


PLACIN 
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SUBJECT: SUBR@UTINE TABL2 


PURPQSE: SUBRGZUTINE TABL2 inputs tables of the form 
F = £(X,¥) square matrix, table format 2Q. 


METH@D: SUBR@UTINE PLACIN is called to read and count the 
Y array, then again to read and count the X array. 
Then for each input Y point PLACIN is called to 
_ read NX points of the dependent array. 


USAGE: CALL TABLE ~2 (T, MAXX, MAXY, NX, NY) 
T -- Table to be read 
MAXX - Maximum number of points in X array 


MAXY - Maximum number of points in Y array 


Output NX - Actual number of X array points input 
Output NY - Actual number of Y array points input. 
ST@RAGE: The table T is stored as follows: 


Address (1) —  Toyayxy  Tmaxxt1)7 7 7 (waxx+Maxy) 


Content Y 


(1) ~ 7 Yonyy 


Address T (MAXX+MAXY+1)~ ~ 7 (MAXX+MAXY+NX) ~7 (12MAXX+MAXWL) 


- « 


F 
Content (1,2) 


Fay aka) 


Address 


T ((NY+1) MAXX+MAXY+1) 7 ((NY+1)MAXX+MAXY+NX) - - 


T ( (MAXY+1) MAXX+MAXY) 


-— = eo ee 2 = eel 


Content 


F (1 ivy) F (nx, NY) 


T can be specified in the calling program by 
either: 


DIMENSI@N T (MAXX+MAXY+ (MAXX) (MAXY)) or 
COMMON/T/TX (MAXX) , TY (MAXY) , TF (MAXX, MAXY) 


4 
. 


| 
Li 


SUBPR@GRAM: 


PLACIN 
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SUBJECT: SUBRGUTINE TABL3 


PURP@SE: SUBR@UTINE TABLE3 inputs tables of the form 
. Fs £(X,Y,2), tabie format 3. 


fe 


METH@D: SUBR@UTINE PLACIN is called to read and count 
the Z array, the Y array and the X array in 
turn. Then is recailed to read a jependent 
array corresponding to the X values for the 
firet Y and 2; then for the second Y, etc. y; 
This is repeated for each Z value in turn. 


-———— ee 8 


USAGE: CALL TABL3 (T, MAXX, MAXY, MAXZ, NX, NY, NZ) 
T - - Table to be read 
| MAXX - Maximum number of points in X array ; 
MAXY - Maximum number of points in Y array 
| MAXZ - Maximum number of points in Z array 
Output NX - Actual number of X array points input 


Output NY - Actual number of Y array points input 
Output NZ - Actual number of Z array points input 


ST@RAGE: The table T is stored as follows: 


Address T(1) T (MAXX) T (MAXX+1) T (MAXX+MAXY ) | | 


Content . X(1) X (NX) Y(1i) Y (NY) 


Address T (MAXX+MAXY+1) T(MAXX+MAXY+MAXZ) T (MAX¥Y+MAXY4MAXZ+1) 


Content Z(1) Z (N2@) F(1,1,1) 


Address 


T (MAXX+MAXY+MAXZ+NX)- - T (2MAXX+MAXY+MAXZ+1) 


Content F (NX,1,1) F(2,1,1) 


NXSMAXX NYSMAXY, NZSMAXZ 


SUBPRZGRAM: 


T can be specified in the calling program by 
either: 


DIMENSIQN T (MAXX+MAXY+MAXZ+ (MAXX) (MAXY! (MAXZ) ) 


or CQMMON/T/TX (MAXX) , TY (MAXY) , TZ (MAXZ, 
TF (MAXX , MAXY , MAXZ) 


Note: dependent array stored in natural 
FORTRAN order. 


PLACIN 
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SUBJECT: 


PURPYSE: 


METH@D: 


USAGE: 


Output 


Output 


ST@RAGF: 


SUBPR@GRAM: 


SUBRGUTINE TABL22 | 


SUBRGUTINE TABL22 inputs tables of the form D0 
F = £(¥,Y) skewed matrix, table format 2K. 


SUBRZUTINE PLACIN is called to read and count 
the Y array. Then for each Y value in turn, 
PLACIN reads and counts and X array and the 


corresponding F values. ‘Th. various X arrays 


need not have the same number of input points. 


CALL TABL22 (T,MAXY,MAXX,NY,NX) 


T = 


MAXY 


NX 


T(1) 
¥(1) 


Table to be read 

Maaimum of points in Y array 

Maximum of points in X array 

Actual number of Y array points input 


Array with actual numbers of X array 
points input 


T (MAXY) T (MAXY+1) T (MAXY+MAXX ) 


Y (NY) X(i,L) % (NX (1) ,1) ; 


{ (MAXY+MAXX11) - - T(MAXY+ (MAXY) (MAXX),1) - - 


X(1,2) X (NX (NY) , NY) 
T (MAXY+ (MAXX) (MAXY)+1 - - T (MAXY+MAXY (2MAXX) ) > 
F(1,1) F (NX (NY) , NY) 


T and NX are svecified in the calling program 
by either: 


DIMENSI@N T(MAXY+MAXY (2MAXX)), NX (MAXY) 


or CQ@MM@ON/T/TY (MAXY), TX (MAXX,MAXY), TF (MAXX, MAXY) 


Note: 


X and Y have been reversed in the actual 
coding of this routine, but to be consistent 
with the write-ups for the other TABL 
routines, the outer loop variable is 

called Y in this write-up. 


PLACIN 
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SUBJECT: 


SeCTION III 


DESCRIPTION OF FUNCTIONS, AND LISTINGS 


FUNCTIQN AREA 


PURP@SE: FUNCTIQ@N AREA computes flow area 
’ RWT 
METHQD : AREA = p-ya- = sq. ft. 
D 
USAGE: A = AREA (W, T, P, Vv, Cy) 
W Inlet mass flow lb/sec 
T Nozzle throat static temperature 
degrees R 
P Nozzle throat static pressure lb/ft? 
VY Nozzle throat velocity ft/sec 
Cy Hiuscharge coefficient = 1.0 
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SUBJECT: 


PURP@SE: 


METH@D: 


USAGE: 


~ Thermodynamic function set; H@FT, T@FH, PROFH, 


C@FH, and C@FHS. 


The thermodynamic function set produces data 
from Keenan and Kaye Gas Tables to the program. 


A polynomial fit of the tabular data is used. 

The appropriate coefficients are det 2rmineéd 

by the checking the range of the input variables. 
H = H@OFT (T,FQA) 

T = T@QFH (H,FQ@A) 

PR = PR@FH(H,FQA) 

H = H@FPR(PR,FQ@A) 

C = C@FH(HT, FPA) 

C = C@FHS (HS,FQA) 

H Enthalpy BTU/1b T temperature degrees R 
FRA fuel/air ratio PR relativepressure 

C sonic velocity ft/sec HT total enthalpy BTU/lb 
HS static enthalpy BTU/lo 

The fuel is assumed to be a hydrocarbon of 


composition, Cy Ho . 
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SUBJECT: 


PURPQSE: 


METHQD: 


USAGE; 


FUNCTIGN K@SUNT . 


FUNCTIQG<N KGUNT counts the number of input 
points on a ten field seven digit input 
card, 


Starting with the right most fieli a test 


is made for a blank field. Upon ~inding 
a non-blank field that fieid number is 
returned as the input count. 
K = K@UNT (CARD, NFTI) 

CARD - Array of input points 


NFI - Maximum fields defined on card 
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SUBJECT: FUNCTI@N DEMAND 


PURP@SE: FUNCTIQN DEMAND computes engine demand, Ad/Acr 
as a function of recovery factor. 


METH@D: DEMAND = Area RFW + Area RFW. +A 
——= -+— = sy 
Cc C An 
where Ay is from Table 7 
BY 
Ac 
USAGE: A@AC = [MAND (RF) 


RF is recovery factor 
Area = AREAF (Mach) 


W = Engine airflow 


W. = Secondary airflow 
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Capture area 
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